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Abstract 

In this paper, we propose two important measures, quantile correlation (QCOR) 
and quantile partial correlation (QPCOR). We then apply them to quantile au- 
toregressive (QAR) models, and introduce two valuable quantities, the quantile 
autocorrelation function (QACF) and the quantile partial autocorrelation function 
(QPACF). This allows us to extend the classical Box- Jenkins approach to quantile 
autoregressive models. Specifically, the QPACF of an observed time series can be 
employed to identify the autoregressive order, while the QACF of residuals obtained 
from the fitted model can be used to assess the model adequacy. We not only demon- 
strate the asymptotic properties of QCOR, QPCOR, QACF, and PQACF, but also 
show the large sample results of the QAR estimates and the quantile version of the 
Ljung-Box test. Simulation studies indicate that the proposed methods perform well 
in finite samples, and an empirical example is presented to illustrate usefulness. 

Keywords and phrases: Autocorrelation function; Box- Jenkins method; Quantile correla- 
tion; Quantile partial correlation; Quantile autoregressive model 
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1 Introduction 



In the last decade, quantile regression has attracted considerable attention. There are 
two major reasons for such popularity. The first is that quantile regression estimation 



( Koenker and Bassett 



19781 ) can be robust to non-Gaussian or heavy-tailed data. In ad- 
dition, it includes the commonly used least absolute deviation (LAD) method as a special 
case. The second is that the quantile regression model allows practitioners to provide more 
easily interpretable regression estimates obtained via various quantiles r e [0,1]. More 
references abou t quantile regre ssion estimations and interpretations can be found in the 



seminal book of 



Koenker! ( 120051 ) . Further extension of quant ile regression to various m odel 



and data st ructures hav e been found in the literature, e.g. 



count data, 



Machado and Silval (I2005h for 



Mu and He] (120071 ) for power tran sformed data, 



Wang and Wangi ( 120091 ) for survival analysis, 



ie and Liand (120001) a nd 



( 120091 ) for regression wit 
augmented factors, and 



i measurement errors, 



Peng and Huanel (120081) and 



Kai et al 



Ando and Tsay 



Wei and Carroll 



(1201 if ) for regression with 



(1201 ll ) for semiparametric varying-coefficient partially 
linear models, among others. 

In addition to the regression context, the quantile technique has been employed to 



the field of time series; see, for example, 



for autoregressive 



( AR) models, 



Koul and Saleh 



Jl995J ) and 



Cai et al. 



(120121 ) 



Ling and McAleerl ( 120041 ) for unstable AR models, and 



Xiao and Koenkerl (120091 ) for g eneralized autoregressive c onditional heteroscedastic (GARCH) 



models. It is noteworthy that 



Koenker and Xiaol (120061 ) established important statistical 



properties for quantile autoregressive (QAR) models, and suggested a modified Bayesian 
information criterion (BIC) to select the order of QAR models. Their findings have ex- 
panded the classical AR model into a new era, which motivates us to extend the classical 
Box- Jenkins' approach (i.e., model identification, model parameter estimation, and model 
diagnostics) from AR to QAR models. In the classical AR model, it is known that model 
identification usually relies on the partial autocorrelation function (PACF) of the observed 
time series, while model diagnosis commonly depends on the autocorrelation function 
(ACF) of m odel residuals. D etailed illustrations of model identification and diagnosis can 



be found in 



Box et al. 



(120081 ). 



The aim of this paper is to introduce two novel measures to examine the linear and par- 
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tial linear relationships between any two random variables for the given quantile r G [0, 1]. 
We name them quantile correlation (QCOR) and quantile partial correlation (QPCOR). 
Based on these two measures, we propose the quantile partial autocorrelation function 
(QPACF) and the quantile autocorrelation function (QACF) to identify the order of the 
QAR model and to assess model adequacy, respectively. It is noteworthy that the appli- 
cation of QCOR and QPCOR is not limited to QAR models. They can be used broadly 
as the classical correlation and partial correlation measures in various contexts. 

The rest of this article is organized as follows. Section 2 introduces QCOR and QP- 
COR. Furthermore, the asymptotic properties of their sample estimators are established. 
Section 3 obtains QPACF and its large sample property for identifying the order of QAR 
model. In addition, the autoregressive parameter estimator and its asymptotic distribution 
are demonstrated. Moreover, QACF and its resulting test statistics, together with their 
asymptotic results, are provided to examine the model adequacy. Section 4 conducts sim- 
ulation experiments to study the finite sample performance of the proposed methods, and 
also presents an empirical example to demonstrate usefulness. Finally, we conclude the 
article with a brief discussion in Section 5. All technical proofs of lemmas and theorems 
are relegated to the Appendix. 



2 Correlations 



2.1 Quantile correlation and quantile partial correlation 

For random variables X and Y, let Q t ,y be the rth unconditional quantile of Y and Q T y(X) 
be the rth quantile of Y conditional on X. One can show that Q T ^(X) is independent 
of X, i.e. Q Ti y(X) = Q t ,y with probability one, if and only if the random variables 



has been used by 



He and Zhu 


(2003 


) and 


Mu and He 


(2007) 



define the quantile covariance given below. For < r < 1, define 

qcov r {F,X} = cov{/(F - Q r , y > 0),X} = E{^ T (Y - Q T>Y ){X - EX)}, 
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where the function ip T {w) = r — I{w < 0). Subsequently, the quantile correlation can be 
defined as follows, 

yorAY.X) - <**"M X > E{4< T (Y - Q T y)(X - EX)} 

Vvar{^ r (y - Q T y)}var(X) \/(t - t>* 

where o\ = var(X). 

In the simple linear regression with the quadratic loss function, there is a nice rela- 
tionship between the slope and correlation. Hence, it is of interest to find a connection 
between the quantile slope and qcov r {Y, X}. To this end, consider a simple quantile linear 
regression, 

(ao, bo) = argmin E[p T (Y — a — bX)], 

a,b 



in wh ich one attempts to approximate Q T y{X) by a linear function ao + b$X (see 



Koenker 



20051 ). where p T {w) = w[r — I{w < 0)]. Then, we obtain the relationship between b and 



qcor r {F, X} given below. 

Lemma 1. Suppose that random variables X and Y have a joint density and EX 2 < oo. 
Then the values of (ao, bo) are unique, and the quantity bo = if and only if the quantile 
correlation qcor T {Y, X} = 0. 

It is noteworthy that the proposed quantile covariance here does not enjoy the symmetry 
property of the classical covariance, i.e., qcov r (Y, X) ^ qcov r (X, Y). This is because 
the first argument of the quantile covariance or the quantile correlation is related to the 
rth quantile, while the second argument is the same as that of the classical covariance. 
Accordingly, qcor r (F, X) ^ qcor T (X, Y). 

Suppose that a quantile linear regression model has the response Y, a q x 1 vector 
of covariate Z, and an additional covariate X. In the classical regression model, one can 
construct the partial correlation to measure the linear relationship between varia bles Y 



and X after adjusting (or controlling) vector Z (e.g., see 



Chatterjee and Hadi . 



20061 ). This 



motivates us to propose the quantile partial correlation function. To this end, let 
{cti,P[) = argmin£(X-a-/3'Z) 2 
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where (a, (3')' is a vector of unknown parameters. Accordingly, a% + (3[Z is the linear effect 
of Z on X. Next, consider 

(a 2 ,(3' 2 ) = argmin£[p r (T -a- f3'Z)}. 

As a result, a 2 +(3' 2 Z is the linear effect of Z on the quantile Y (i.e., the linear approximation 
of Q TiY (Z)). It can also be shown that E{X-a 1 ~P[Z) = 0, E[ifj T (Y - a 2 - (3' 3 Z)} = and 
E[ip T (Y — a 2 — $$Z)Z] = if the random vector (X, Y, Z')' satisfies the conditions stated 
in the forthcoming Lemma [2j Using these facts, we define the quantile partial correlation 
as follows, 

cov{i> T (Y -a 2 - /3' 2 Z),X - ai - P[Z} 



qpcor r {F,A|Z} 



^/var{^ T (F - a 2 - /^Z)}var{X - a x - f3[Z} 
E[ip T (Y - a 2 - (3' 2 Z)(X -a,- (3[Z)] 

y/{T-T*)E{X- ai -P[Z)* 

E[jj T ( Y-a 2 -(3' 2 Z)X} 

r - r2 )4| Z 



(2.2) 



where o~ 2 x ^ z = E(X — ct\ — f3[Z) 2 . This indicates that the covariate X has no additional 
linear contribution to the quantile response Y iia 2 +f3' 2 Z = a^+^Z+^X with probability 
one, where 

(« 3 , £3,73) = argmin£[p T (r - a - (3'Z - jX)}. 
This leads to the following lemma. 

Lemma 2. Suppose that the random vector (X, Y, Z')' has a joint density with EX 2 < 00 
and E\\Z\\ 2 < 00, where \\ ■ || is the Euclid norm. Then (a3,/3 3 ,73) = (a 2 , (3' 2 ,0) if and only 
if the quantile partial correlation qpcor T {Y, X\Z} = 0. 

Since the true qcor r and qpcor r are often unknown in practice, we introduce their sample 
versions given below. 

2.2 Sample quantile correlation and sample quantile partial cor- 
relation 

Suppose that the data {{Yi, X iy Z^)', i = l,...,n} are identically and independently gener- 
ated from a distribution of (Y, X,Z')' . Let Q t ,y = inf {?/ : F n (y) > r} be the sample rth 



quantile of Y ly Y m where F n (y) = n _1 EILi H^i < y) is the empirical distribution func- 
tion. Based on equation ( 12. ip . the sample estimate of the quantile correlation qcor r {F, X} 
is defined as 



1 1 ~ 

qcor T {y, X} = ■ - V MYi ~ Q.,y)(X - X), (2.3) 

where X = n" 1 X, and a| = n" 1 £" = i(** - ^) 2 - 

To study the asymptotic property of qcor r {Y, X}, denote fy(-) and /y|x( - ) as the 
density of Y and the conditional density of Y given X, respectively. In addition, let 
fix = E(X), fi x \Y = E[f Y] x(Qr,Y)X]/f Y (Q TiY ), E n = E(X - fi x f - a x , 

E 12 = - Q t ,y)(X - fx x \ Y )} 2 - [qcov r {F,X}] 2 , 

E 13 = E[ij T (Y - Q T , Y )(X - fi x \Y){X - fix) 2 ] - 4 • qcov r {F,X}, 

and 



fix 



T — T 2 



Sntqcov^X}) 2 S 13 ■ qcov r {F, X} S 12 
4a 6 v ai + a 2 



'X u x u x 

where o~ x is defined as in the previous subsection. Then, we obtain the following result. 

Theorem 1. Suppose that E(X) 4 < oo and there exists a n > such that the density 
jfy(-) is continuous and the conditional density /y|x( - ) is uniformly integrable on [Q Tj y — 
ti", Qt,y + 7r] . T/ien 

(?corT T {y, X} - gcorr r {F, X}) -> d X(0, fi^. 

To apply the above theorem, one needs to estimate the asymptotic variance Q±. To 
this end, we employ a nonparametric approach, such as the Nadaraya- Watson regression, 
to estimate the function m(y) = E(X\Y = y), and denote it as fh(y). We further assume 
that the random vector (X, Y) has a joint density, which leads to fix\Y = E(X\Y = Q t .y)- 
Accordingly, we obtain the estimate, fix\Y — "^(Qt,y)j where Q t ,y is the rth sample 
quantile of {Y ly ...,Y n }. Finally, the rest of quantities, fix, & x , qcov T {Y, X}, En, E 12 , 
and E13 contained in fl± can be, respectively, estimated by 'fix = X = n~ l Y^i=\ X%, & x = 
n' 1 ELi^-^x) 2 , qcov T {Y, X} = n~ l Y!U ^-Q^XX-X), E n = £™ =1 (X- 
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fix) 4 - a x , E 12 = n-^tAMYi ~ Qr,Y){Xi - fi x \ Y )} 2 - [q5ov r {Y,X}] 2 , and E 13 = 
n _1 Y^=i ^riXi — Qr,Y)(Xi — fix\Y)(Xi — fix) 2 — o\ qcov r { Y, X } . As a result, we obtain 
the estimate of fii, and denote it by fii. 

We next estimate the quantile partial correlation qpcor r {Y, X}. Let 

n n 

(Si, = argmin — a — /3'Zj) 2 and (a 2 , /3 2 ) = argmin p T (Y; — a — /3'Z,j). 

Based on equation (12. 2p . the sample quantile partial correlation is defined as 

1 1 n 

qp^? r {F,X|Z} = - • -y>r(*; - «2 - /3 2 Z,)X,, (2.4) 

V / (- - - 2 )^,z n ^ 

where a| |z = n~ l £™ =1 (Xi - «i - ^Z 4 ) 2 . 

To investigate the asymptotic property of qpcor r {Y, -AT|Z}, denote the conditional 
density of Y given Z and the conditional density of Y given Z and X by /y|z( - ) and 
/r|z,x(-); respectively. In addition, let 6 1 = (c*i, /?()'> #2 = (a2,/? 2 )', z * = (1» z 0'> S 2i = 
£'[,/V|z,x(# 2 Z*)XZ*], ^22 = E[f Y \z(8 2 Z*)Z*Z*'], E 20 = E 21 E 22 , E 23 = E(X— 9[Z*) 4 — a x \ z , 

E 24 = E[MY - 9 2 Z*)(X - E 20 Z*)] 2 - {E[ifj T (Y - 9 , 2 Z*)X}} 2 1 

E 25 = E[,f> T {Y - 8 2 Z*)[X - E 20 Z*)(X - ^Z*) 2 ] - a 2 ■ E[^ T {Y - 9' 2 Z*)X], 



and 



SI 



2 



r — t 2 



£ 23 (£tyv(y - # 2 Z*)X]) 2 E 25 ■ E[^ T (Y - 9' 2 Z*)X] E 



+ 



21 



4rr 6 rr 4 rr 2 

^°X|Z °X|Z °X|Z 



where cci, 0i, a 2 , /3 2 and crL z are defined as in the previous subsection. Then, we have 
the following result. 

Theorem 2. Suppose that E 21 < 00, < E 22 < 00, EX 4 < 00, E\\Z\\ 4 < 00, E(Z*Z*') > 
0, and there exists a tt > such that fy\z{9' 2 Z* + ■) and fY\z,x(9' 2 Z* + •) are uniformly 
integrable on [—7T,7i]. Then 



n[qp7b~r T {Y,X\Z} - qpcor T {Y, X\Z}} ^ d JV(0,fi 2 ). 

To estimate the asymptotic variance Q 2 given in Theorem 2, we consider Y* = Y — 6' 2 Z* 
and qcov r {y*,X} = E[ip T (Y — 8' 2 Z*)X]. In addition, assume that the random vector 
(Y,X,Z'Y has a joint density. We then have that E 2 i = E[f Y * \z,x{0)XZ*} = f Y *(0) ■ 



E[XZ*\Y* = 0], £ 22 = f Y *(0)-E[Z*Z*'\Y* = 0], and £ 20 = E[XZ*'\Y* = 0}{E[Z*Z*'\Y* = 
0]} -1 , where /y*(-) is the density of Y*. Applying the same nonparametric technique as 
that used for estimating hx\y in Theorem 1, we could estimate each of the vector and matrix 
components in va.i(y) = E[XZ*\Y* = y] and m 2 (?/) = E[Z*Z*'\Y* = y], respectively, from 
the data {(Y*, X h Z'J = (Y i -6 2 Z*,X h Z' i ),i = l,...,n}, where 6 2 = (a 2 ,3 2 ) / . Accordingly, 
S 20 = S^S^ 1 = m' 1 (0)[m 2 (0)]~ 1 . Subsequently, the rest of quantities involved in f2 2 , o" 2 ^, 
qcov r {F*,X}, S 23 , S 24 , and E 25 can be, respectively, estimated by cr^| Z = n~ l Y^i=ii-^i ~ 
a 1 -^[Z l ) 2 , qcov T {Y*, X} = n' 1 J™ =1 ^(^-^Z*)^, S 23 = n' 1 Er=i(^-#i Z I) 4 -^|z> 
S 24 = n~ l Er=i[^(^-^Z*)(X,-S 20 Z*)] 2 -[q56V T {r*, X}] 2 , and S 25 = n~ l £™ = i M Y i~ 
6 2 Z*)(Xi-Tj 2 oZ*)(Xi-6[Z*) 2 -a^^ z -c^ov T {Y* , X}. Consequently, we obtain the estimate 
of f2 2 , and denote it by f2 2 . 

It is noteworthy that the quantile correlation and quantile partial correlation can be 
broadly used as the classical correlation and partial correlation in regression analysis (e.g., 
variable selections), although our focus is on quantile autoregressive models. 



3 Quantile autoregressive modeling 

Suppose that {y t } is a strictly stationar y and ergodic time s eries, and Tt is the a-field gen- 



erated by {y t ,yt-i, ■■■}■ We then follow 



Koenker and Xiaof s (2006) approach and present 



QAR models; i.e., conditional on J-~t-i, the rth quantile of y t has the form of 

Q T {yt\Ft-i) = <f>o(r) + 0i(r)y t _i H h (f)p(r)yt-p for < r < 1, (3.1) 

where <pi{-)s are unknown functions mapping from [0, 1] — > R. Following the Box- Jenkins' 
classical approach, we next introduce the QPACF of a time series to identify the order of 
a QAR model, and then propose using the QACF of residuals to assess the adequacy of 
the fitted model. 



3.1 Model identification and estimation 

For the positive integer k, let z tife _i = (y t -\, Vt-k+i)', («i, P{) = argmin Q /3 E(y t _ h - a - 
P' z t,k-i) 2 , and (a 2 , f3' 2 ) = argmin Q /3 E[p T (y t — a— f3'z ttk _i)], where the notations (ai, (3[) and 
(a 2 , (3 2 ) are a slight abuse since they have been used to denote the regression parameters 
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in Section 2. From equation ( 12. 2p . we obtain the quantile partial correlation between y t 
and y t -k after adjusting the linear effect Zt,k-ii 

, f | ! E[ip T (y t - a 2 - /02 z *,fe-i)2/*-fc] 

0fcfc,r = qpcor T {y t ,y t -fc|zt,fc-i} = =^7 fl/ Ff » 

^/(t - T 2 )E(y t _ k -a 1 - /3{z t)fc _i) 2 

and it is independent of the time index t due to the strict s tationarity of {yt}- Analogous 



to the definition of the classical PACF (IFan and Yao . 



2003 



Chapter 2), we name <fikk,T to 



be the QPACF of time series {yt}- It is also noteworthy that </>n, T = q c or T {y t , yt-i}- We 
next show the cut-off property of QPACF. 

Lemma 3. If ' <t> p {j) 7^ with p > 0, Ey\ < 00 and E[y t — E (y t \J-'t-i)] 2 > 0, then 4> PP)T 7^ 0, 
and 4>kk,r = for k > p. 

The above lemma indicates that the proposed QPACF plays the same role as that of PACF 
in the classical AR model identification. 

In practice, one needs the sample estimate of QPACF. To this end, let 

n n 

(a 1 ,(3[) = argmin ^ {yt-k-a-ft^k-i) 2 \ {^2,^2) = argmin ^ pAvt-a- (3'zt,k-i), 
a >P t=k+i a ^ t=k+i 

and = n~ l XltLfc+i G/t-fc — ot\ — 0[zt,k-i) 2 - According to ( 12 .4p . we obtain the estimation 

for <f> kk<T , 



'kk,T 



T - r2 ) CT y|z n t=k+l 



1 - 

- Y] iPAyt - a 2 - f3' 2 z t,k-i)yt-k, 



and we term it the sample QPACF of the time series. 

To study the asymptotic prope rty of 4>kk,r, we introduce the following assumption, 



which is similar to Condition A. 3 in 



Koenker and Xiaol (120061 ) . 



Assumption 1. Eyf < 00, E[y t — E(y t \ Ft^i)} 2 > 0, and there exists a n > such that 
/t-i(') i s uniformly integrable on [— tt,tt]. 

Furthermore, let 

e t , T = yt- M r ) - 4>i( T )yt-i M T )yt- P - (3.2) 

By (13. ip . the random variable I{&t >T > 0) is independent of y t -k for any k > 0, and 
(0:2,^2) = (4>o(t), 0i(t), p (r), 0, 0) for k > p. Let ft-i(-) be the conditional density 
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of e t>T on the cr-field F t -x, and z* fc _ 1 = (1, zJ^J' = (l,y t -i, ...,y t -k+x)' ■ Moreover, let A = 

E[yt-kK,k-i\, M = E[f t ^)yt-kzik~il s 30 = £K*_i<* -J. s si = £[/ t -i (oKv^-il, 

and 

^ E(y t 2 ) - gA^Ao + A^g^gA, ^ 
Then, we obtain the asymptotic result given below. 

Theorem 3. For k > p, if Ax < oo, < Ift < oo and Assumption^ is satisfied, then 
4>kk,r = and 

Vn<f>kk,T -*d N(Q,tt 3 ). 

To estimate f2 3 in the above theorem, we first apply the 
method to obtain the estimation of ft-x(0) given below. 

2h 



Hendricks and Koenkerl (Il99ll ) 



f t ^(0) = ^ 

Q T+h {yt\jF t -x) - Qr-hiy^Ft-x) 

where Q T (yt\Ft-x) — 0o( r ) + ft( r )?/t-i + ■ • • + <Pk(T)yt-k is the estim ated rth quantile o: 
lit an d h is the bandwidth selected via appropriate methods (e.g., see 



Koenker and Xiao 



20061 ). Afterwards, we can use the sample averaging to approximate Aq, Ax, £30, £31, 
E(y 2 ), and E(y t ~ k - a x - ftz ijfc _i) 2 by replacing their f t -i{-), a ± , and ft, respectively, 
with /t-i(0), &x and ft. Accordingly, we obtain an estimate of f2 3 , and denote it as Q3. 



In sum, we are able to use the threshold values ±1.96y f^/n to check the significance of 

4>kk,T- 

To demonstrate how to use the above theorem to identify the order of a QAR model, we 
generate the observations yx, 2/200 from y t = + a(ut)yt-x, where <3> is the standard 

normal cumulative distribution function, a(x) = max{0.8 — 1.6%, 0}, and {u t } is an i.i.d 
sequence with uniform distribution on [0, 1]. We attempt to fit the QAR model ( 13. ip with 
r = 0.2, 0.4, 0.6, and 0.8, respectively, to the observed data {yt}- Figure [1] presents the 



sample QPACF 4> k k,T for each r with the reference lines ±1.96y f2 3 /n. We may conclude 
that the order p is 1 when r = 0.2 and 0.4, while p is when r = 0.6 and 0.8. 

After the order p of model (13.1 ft is correctly identified, we subsequently fit the selected 
model to data. Let <fi = (</>o, ft, 4> p )' be an any parameter vector in model (13.11) and 
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4>( T ) — ((0o( r ); 0i( r )? •••?0p( r )) / be the true value of <f>. Consider 

n 

0(r) = argmin p T (y t - (f>'z* t p ), 

4> t=p+i 

where z* p = (l,z' tp )' = (1, y t _ u y t - P )' ■ In addition, let E 40 = £[z* p z*' p ], E 41 = 
-E , [/t-i(0)zJ , p z^' ;p ], and f2 4 = (r — r 2 )E 4 " 1 1 E4 E 4 " 1 1 . We then obtain the following asymp- 
totic property of the estimated parameter. 

Theorem 4. If < E 4 i < oo and Assumption^ is satisfied, then 
v^{0(r) - 0(r)} ^ d iV(0,fi 4 ). 



Koenker and Xiaol (120061 ). although we 



The above result is similar to that of Theorem 2 in 
make different assumptions. The Jl 4 in the above theorem can be estimated by applying 
the same techniques used for the estimation of Q3. 

3.2 Model diagnostics 

For the errors {e tjT } defined in (13. 2p . we employ equation (12. ip and the fact that Q T ,e t , T = 0, 
and obtain QACF between {e ttT } and {e t -k, T } as follows, 

E{tp T (e t)T )[e t -k,T ~ E(e t , T )]} 



Pk, 



y/(r - r 2 )a 2 e 



where a\ = var(e fiT ). Suppose that the QAR model is correctly specified. We can show 
that pi~ T = for k > 0. Hence, we are able to use pk )T to assess the model fit. In the 
sample version, we consider the residuals of the QAR model, 



e t , T = Vt~ M T ) ~ M T W-i 9p{VVt- 



for t = p + 1, ...,n, and e^ r = for t — 1, ...,p. It can be verified that the rth empirical 
quantile of {e ttT } is zero. Based on this fact and equation (12. 3p . we obtain the estimation 

of Pk,T, 

1 1 n 

rk,r = 7= == • - Y] 1p T (e t ,r)(et-k,r - He), 

y/(r-r 2 )a 2 e n t=fc+i 

where k is a positive integer, Jl e = n~ l X)t==fc+i ^*,tj = n ~ 1 Y^t=k+i(^i>r ~ Pe) 2 , and the 
rth empirical quantile of {et, T } is zero. We name Tk, T the sample QACF of residuals. 
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Adapting the classical linear time series approach (ILil . 120041 ). we examine the signif- 
icance of {rfc jT } individually and jointly. For the given positive integer K, let e,t-i,K — 
(e t -i, T , ■■■,e t -K,r)', S 50 = E[et-i,Kzt'p], S 51 = E[/ t _i(0)e t _i i KZf ' p ], and 

^5 = - ^{-E'( e t-i,K e t-i,i<") + ^si^i^o^/E^ — £ 5 iE 41 1 Eg — E 50 X! 41 1 £g 1 }. 

Then, we obtain the asymptotic distribution of R T = (r l r , ...,r K)T )' given below. 

Theorem 5. Assume that < X 41 < oo, £ 51 < oo, and Assumption^ holds. We then 
have 

Applying the same techniques as used in the estimate of f2 3 , we are able to estimate the 
asymptotic variance and denote it Q5. In addition, let the k-ih diagonal element of Q5 



be f2 5 fc. Then, one can employ r^^j y VL^ to examine the significance of the fc-th lag in 
the residual series. 

To check the significance of R T jointly, it is natural to consider the test statistic 
R' T Vt§ l R T . However, may not be invertible. Hence, we approximate ^5 by Ik — 
cr~ 2 E 50 Sj 1 ^'50' w hich holds under the assumption that {e tjT } is an independent and iden- 
tically distributed (i.i.d.) sequence and /(_i(0) is a constant. The resultin g matrix is idem- 



poten t and has rank K — p. This allows us to obtain a Box-Pierce type (IBox and Pierce 



19701 ) test statistic, 



Q BP (K)=nJ2 

5=1 



K 

2 



r ■ 

],T1 



which follows an approximately chi-squared distribution with K — p degrees of freedom, 
xk- p - Accordingly, Qbp(K) can be used to test the significance of pi >T to pk, t jointly. 



4 Simulations and an empirical example 
4.1 Simulation studies 

We conduct five simulation experiments to assess the finite-sample performance of the 
proposed methods. Specifically, the first simulation experiment is for the sample quantile 
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correlation and the sample quantile partial correlation proposed in Section 2, and the last 
four experiments are, respectively, for identification, estimation, and diagnosis as intro- 
duced in Section 3. In all experiments, we conduct 1000 realizations for each combination 
of sample sizes n = 50, 100, and 200 and quantiles, r = 0.25, 0.50, and 0.75. 

In the first simulation experiment, we generate the i.i.d. samples {(Xj,Yj, Zj),i = 
1, ...,n} from the following multivariate normal distribution, 

/ 1.0 0.5 0.5 \ ' 



(X, Y,Z)~N< 



0. 



0.5 1.0 0.5 
\ 0.5 0.5 1.0 / t 
After algebraic simplification, we obtain that 

qcor T {F,X} = 0.5 exp{-0.5 (r)] 2 }/ y/(r - t 2 )2tt, 

and qpcor T {Y, X\Z} = qcov T {Y,X}/\/3, where $(•) is the cumulative standard normal 
distribution. Tables [T] and [2] present the bias (BIAS) and estimated standard deviation 
(ESD), respectively, of the sample quantile correlations qcor T {Y, X} and the sample quan- 
tile partial correlations qpco? T {Y, X\Z}, calculated from 1000 realizations. 

To estimate the asymptotic variances Qx and Q2, we mainly need to estimate the 
quantities hx\y and £ 2 o, addressed in Subsection 2.2. To this end, we emplo y the Nadaraya - 



Wat son approach with the tw o bandwidth selection methods proposed by 



and 



Bofingerl (119751 ) 



Hall and Sheatherl (119881 ). respectively, which are given below. 



Kb = n 



1/5 f 4.50 4 ($- 1 (r)) 



1 [2($-W + l] 2 



1/5 



and h HS 



n -l/3^2/3 



1.5(f) 2 (^(t)) 



1/3 



where </>(•) is the standard normal density function, z a = $ _1 (1 — a/2), for the construction 
of 1 — a confidence intervals, and a is se t to 0.05. Furthermore, we consider two more 



bandwidths, 0.6/is and 3hns, suggested by 



Koenker and Xiaol (120061 ). In sum, we have four 



bandwidth choices. The resulting asymptotic variance estimates, Qi and yield their 
correspondingly asymptotic standard deviations (ASDs) given in Tables [T] and [2J Both 
tables indicate that the ASDs are close to their corresponding ESDs even when n = 50, 
and they become smaller as the same size gets larger. In addition, biases are close to zero, 
and decrease as the sample size increases. Moreover, all four bandwidths lead to similar 
results, although 3hns is slightly better than the others. 
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The second Monte Carlo experiment studies the performance of QPACF for identifying 
the order of the QAR models, we generate the data from the following process, 

yt = 0.1 + 0.5^_i + et, (4.1) 

where {e t } is an i.i.d sequence with standard normal distribution. Under the above setting. 



i t can be shown that = 1. We then employ the approach of 



Hendricks and Koenker 



(jl99ll ) with the four bandwidths used in the first experiment to estimate the density 
function, / t _ 1 (0). This allows us to further estimate the variance matrix 3 in Theorem [21 
(see Subsection 3.1). Table [3] presents the bias and estimated standard deviation of 4>kk,r 
at k — 2, 4, and 6. It shows that biases are small even when n = 50, and the ESDs are 
close to the ASDs as well as their theoretical value 1 / y/n. 

The third simulation experiment investigates the finite-sample performance of the QAR 
estimates. We use the same data generated from (14.11) . and then fit it with the QAR 
model (13.11) with p — 1. In addition, we employ the same approach as given in the second 
experiment to estimate f t _i(0). As a result, the variance matrix f2 4 in Theorem H] can be 
estimated (see Subsection 3.1). Table H] presents the biases, estimated standard deviations, 
and asymptotic standard deviations of parameter estimates </>o(t) and <Pi(t). It shows that 
biases are close to zero even when the sample size is as small as n = 50. In addition, 
the ESDs are close to the ASDs, and both of them decrease as the sample size increases. 
Moreover, there is no discernible difference among the four bandwidths, although 3hns 
often yields the smallest ASD. 

The fourth simulation experiment examines the finite-sample performance of the sample 
QACF of residuals individually via the asymptotic result in Theorem [5] All settings are 
the same as those in the third experiment. Table [5] presents the biases, estimated standard 
deviations, and asymptotic standard deviations of r^^ at k — 1, 3, and 5. Apparently, 
biases are small and the ASDs are close to their corresponding ESDs. 

Finally, the fifth experiment studies the approximate test statistic Qbp{K)- To this 
end, we generate data from the following process, 

yt = 0.5y t _i + <f)y t -2 + e t , 

where {e t } are i.i.d. standard normal random variables. In addition, = corresponds to 
the null hypothesis, while > is associated with the alternative hypothesis. Moreover, 
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the nominal level is 5%. Table E] reports sizes and powers of Qbp{K) with K = 6. It shows 
that Qbp{K) controls the size well, and its power increases quickly when the sample size 
or <p becomes larger. Consequently, the above six simulation studies perform satisfactorily 
and support our theoretical findings. 



4.2 Nasdaq Composite 

This example considers the log return (as a percentage) of the daily closing price on 
the Nasdaq Composite from January 1, 2002 to December 31, 2007. There are 1,235 
observations in total, and Figure [2] depicts the time series plot and the classical sample 
ACF. It is not surprising to conclude that these returns (i.e., log returns) are unc orrected 



and ca n be treated as an evidence in support of the fair market theory. However, 



Veronesi 



( 119991 ) found that the stock markets u nder-react to good news in bad times and over-react 



to bad news in good times. Hence, 



Baur et al 



(120121 ) proposed aligning a good (bad) 



state with upper (lower) quantiles by fitting their stock returns data with the QAR(l) 
type models. This motivates us to employ the general QAR model with our proposed 
methods to explore the dependence pattern of stock returns at a lower quantile (r = 0.2), 
the median (r = 0.5), and an upper quantile (r = 0.8). 

We first fit the returns at the lower quantile (r = 0.2), and then present its sample 
QPACF in Panel A of Figure [3j It shows that lags 1, 2, and 13 are significant, which 
suggests QAR(13) could be considered for model fitting. We then refine the model via the 
backward variable selection procedure at the 5% significance level. The resulting model is 

QoM^t-l) = -0.4114 .0269 + 0.1117o.0482?/t-l + 0.0951o.047l2/t-2 + 0.0992o.o4572/t-13, (4.2) 

where the subscripts of parameter estimates are their associated standard errors. Accord- 
ingly, the above coefficients are all significant at the 5% significance level. In addition, 
the second graph in Panel A presents the sample QPACF of residuals, and no lags stand 
out. This, together with the p-value of Qsp(18) being 0.742, implies that this model is 
adequate. 

We next consider the scenario with r = 0.5. The sample QPACF in Panel B indicates 
that all lags are insignificant. Hence, we fit the following model, 

QoMFt-i) = 0.0036 .oi70- (4.3) 
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The above coefficient is not only small, but also not significant. In addition, none of the 
lags in the sample QACF of residuals in Panel B show significance. Moreover, the p- value 
of <3bp(18) is 0.566. Consequently, the above model is appropriate. 

Finally, we study the upper quantile scenario with r = 0.8. The sample QPACF in 
Panel C exhibits that lags 1, 2, 7, 10 and 15 are significant, and suggests QAR(15) could be 
considered for model fitting. After refining the model via the backward variable selection 
procedure, we obtain 

QosiUtl^t-x) = 0.3988 .o230 - 0.1076 o . O3 7i2/t-i - 0.0825 .o4i22/t-2 

- 0.0790 .o36i2/t-io - 0.0802 .o302^i5, (4.4) 

where all coefficients are significant at the 5% significance level. In addition, the sample 
QACF of residuals in Panel C displays that all lags are insignificant. This, in conjunction 
with the p-value of Qbp(18) being 0.215, indicates that the above model fits the data 
reasonably well. 

Based on the three fitted QAR models, (14. 2p . (14.3)) . and (14. 4p . we obtain the following 
conclusions, (i.) The lag coefficients at the lower quantile (r = 0.2) are all positive. 
This indicates that if the returns in past days have been positive (negative), then today's 
negative return is alleviated (even lower). It also implies that stock markets under- react 
to good news in bad times, (ii.) The lag coefficients at the upper quantile (r = 0.8) are all 
negative. This shows that if the returns in past days have been negative (positive), then 
today's positive return (r = 0.8) is even higher (dampened). As a result, stock markets 
over-react to bad news in good times, (iii.) The intercept at (r = 0.5) has a small value 
and is insignificant at the 5% significance level. Thus, the conditional median of returns 
is almost zero as we expected. In addition, equation (14. 3 p indicates that today's return is 
not affected by the returns of recent past days. Although we only report the results of the 
lower and higher quantiles at r = 0.2 and r = 0.8, our studies yield the same conclusions 
across various lower and upper quantiles. In sum, our proposed methods support Veronesi's 



( 119991 1 equilibrium explanation for stock market reactions. 
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5 Discussion 



In quantile regression models, we propose the quantile correlation and quantile partial 
correlation. Then, we apply them to quantile autoregressive models, which yields the 
quantile autocorrelation and quantile partial autocorrelation. In practice, the response 
time series may depend on exogenous variables. Hence, it is of interest to extend those 
correlation measures to the quantile autoregressive model with the exogenous variables 
given below. 

p <? 

i=l j=l 

where is a vecto r of time series, and 4>i{r) and /3j(r) are functions [0,1] — > R, see 



Galvao et al. 



(120121 ) . In addition, the application of the proposed correlations to the quan- 
tile regression model with autoregressive errors is worth further investigation. Clearly, 
the contribution of the proposed measures is not limited to those two mod els. For ex- 



Fan and Lv 



2008 



Wang 



20091 ) in quantile 



ample, variable screening and selection (e.g. 
regressions are other important topics for future research. In sum, this paper introduces 
valuable measures to broaden and facilitate the use of quantile models. 



Appendix: technical proofs 

Proof of LemmaUl For a,b 6 R, denote the function h(a,b) = E[p T (Y*)}, where Y* = 
Y — a — bX. We first show that h(a, b) is a continuously different iable function and has 
derivatives, 

d J^ = -ElMY*)] = P(Y*<0)-T and *M = -E[MY')X]. 
For u ^ 0, 



p T (u - V) - p T {u) = —V1p T (u) + 



/ [I(u < s) - I(u < 0)]ds 
Jo 



-vip T (u) + (u- v)[I(0 > u > v) - 1(0 < u < v)\, (A.l) 
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see iKoenker and Xiaol (120061 ). This, together with Holder's inequality and the fact that 



|y*|/|X| is a continuous random variable, leads to 

\^[h(a, b + c)- h(a, b)] + E[tp T (Y*)X]\ 

= \-E[p T (Y* - cX) - p T (Y*)} + E[iP T (Y*)X]\ 
c 

= \-E{(Y* - cX)\I(0 > Y* > cX) - 1(0 <Y* < cX)}}\ 
c 

< E[\X\I(\Y*\ < \c\ ■ \X\)] < (EX 2 ) 1/2 [P(\Y*\/\X\ < \c\)Y/\ 

which tends to zero as c — > 0. Accordingly, dh(a,b) jdb is obtained. Analogously, we have 
dh(a,b)/da. By Holder's inequality, we can further prove the continuity of dh(a,b)/db. 
Moreover, the continuity of both X and Y implies that dh(a, b)/da is a continuous function. 
It is noteworthy that h(a, b) is a convex function with lim (1 2 +fe 2_ i . 00 h(a, b) = +oo. This, in 
conjunction with the above results, demonstrates that the values of a® and b satisfy 

E[tfj T (Y - a - b X)} = and E[ip T (Y - a - b X)X] = 0. (A.2) 

We next show the uniqueness of (ao,&o)- Suppose that there is another pair of values 
(di, &i) such that h(ai, bi) = h(a ,b ) = argmin a b E[p T {Y — a — bX)}. Let Y = Y — a — b X 
and ^ = (ai - a ) + (&i - b )X. Then, by flATjl and f TA~2|) . 



= h{a u h) - h{a , b ) = E[p T (Y - - p T (Y )] 
= -E[& T (Y )] + E[(Y - 0/(0 >Y > 0] + E[(£ - Y o )I{0 <Y < 0] 
= E[(Y - 0/(0 >Y > 0] + E[(£ - F )/(0 <Y < 0]. 

Note that both (Y - 0/(0 > F > and (f ~ Y )I(0 < F < are nonnegative 
random variables, and Yq — £ is a continuous random variable. Thus, with probability one, 
1(0 > Y > = /(0 < Y < = 0, which implies (ai, 6i) = (a , & )- 

Finally, if b = 0, then (1A.2j) leads to a = Q T y and qcov T {y, X} = E[ip T (Y — ao — 
b X)X] = 0. On the other hand, if qcov T {F, X} = 0, then equation (1A.2j) with (a ,b ) = 
(Qr,y,0) holds. By the uniqueness property, bo = 0, which completes the proof. 

Proof of Lemma\B Let Y* = Y - a 2 - (3' 2 Z and 

(04,^4,74) = argmin E[p T (Y* - a - (3'Z - 7X)]. 

c*,/3,7 
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Since the random vector (X, Y*,Z')' has a joint density, we apply similar techniques to 
those in the proof of Lemma [T] to show that 

E[if> T (Y*)] = 0, E\4, T {Y*)Z] = 0, (A.3) 

and the values of a 4 , {3' 4 and 74 are unique and satisfy 

E[^ T {Y* - « 4 - P' 4 Z - 74 X)(1, Z',X)'} = 0, (A.4) 

where is a (q + 2) x 1 zero vector. 

From ([Q)l . if (a 4 ,/3 4 ,7 4 )' = 0, then qcov r {Y"*,X} = £[^ T (T*)X] = 0. On the 
other hand, qcov r {F*,X} = 0, together with (IA.3[) . implies that equation (IA.4I) with 
(a4, /3 4 , 74)' = holds. Accordingly, we have shown that qcov r {y*, X} = if and only if 
(a 4 ,/3 4 ,7 4 )' = 0. Based on the definitions of (a 2 , /3 2 ) and (a3,/?3,73) in Subsection 2.1, we 
further have that a 4 = a 3 — a 2 , /3 4 = /3 3 — /3 2 , and 74 = 73. Finally, using the fact that 
qpcor r {Y, X|Z} = qcov r {F*,X}/ — t 2 )^^ completes the proof. 

Proof of Lemma For k = p, let 

(a3,/? 3) T3) = argmin£[p r (^ - a - /?'z tiP _i -7^)]. 

It is noteworthy that (03, ,$3,73) = (0o(t), 0i(t), <p p (r)). Since p (r) 7^ 0, we apply 
Lemma [2] and are able to show that <p pp , T 7^ 0. 

Let e t>T = y t - (j> (r) - (j)x(r)yt-i (j} p {T)y t - p . By (JSH]), I(e t)T > 0) is independent 

of y t -k for any k > 0. In addition, (a 2 , /3 2 ) = (0o(t), •••> <Ap( r ); 0') f° r k > p, where 

is (k — p) x 1 vector. Hence, 0fcfc jT = for > p. 

Proof of TheoremUi For m ^ 0, we have that 

I(u - v < 0) - I(u < 0) = JO > m > 0) - I(v <u<0). 

Using this result, we then obtain 

-1 n 1 n 1 1 n 

-y>rP5-§T,y)(*i-*) = - V^(^-g T , y )X,+-A n -X- y>r(^-§T,y), (A.5) 

i=l i=l i=l 

where A n = YT i= i 9r(Yh Qt,y, Qr,v) X i and 

9r(Xi, Qt,Y, Qt,y) 

= MY - Q t ,y) ~ MY - Qt,y) = -[I(Yi < Q T , Y ) - I(Y < Q tX )} 

= I(Qt,Y ~ Qr,Y <Y- Q T y < 0) - I(Q T ,y - Q T ,Y >Y~ Q T y > 0). 
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It can be shown that 



1 _ 1 _ f 1 1 

|-y>rtti - Sr,y)| = |r - - E J (^ - §r,y)l = |t - ^| < 
n n n n 

i=l i=l 

This, together with the law of large numbers, implies the last term of (1A.5|) satisfying 



1 n 

X--J2 M Y i ~ Qr,r) = p (n~ l ). (A.6) 



n 

i=l 



We next consider the second term on the right-hand side of (IA.5I) . For any v G R, 
denote 



where 



1 n 

-= Y^9r{Yi, Q T>Y , Q T , Y + n^v) - E[g T (Y h Q T>Y , Q rX + rC^v^X^Xi 

v i=l 



rQ-r,Y+n-y 2 v 

E[g r {Y h Q T , Y ,Q T , Y + n- 1 / 2 v)\X l ] = - / f Yi \xMdy 



Qt,Y 

and fy^Xii') is the conditional density of Y~j given Xj. Then, by Holder's inequality, we 
have that 

E[Uv)? = E[g T (Y il Q T ,Yi Qt,y + n- l ' 2 v)X l ] 2 

< [PQYt - Q TtY \ < n-^vf^EXt] 1 ' 2 = o(l). (A.7) 

After algebraic simplification, we further obtain 

SUp |£ n (Vi) - £ n (v)\ 
|t>i— v\<5 

1 " 

< sup ^Y,\^M)-9r{v)}X l \+E[\{g T {v l )-g T {v)}X t \\X l ] 

1 U 

= -= IKK) - HrWm + E[\{g T (vl) - grWXiWXl 
where v* takes the value of v + 8 or v — 6. Hence, 

E SUp \£ n (yx) - £ n (v)\ 
\vi—v\<5 

<2 v ^E\{g T (v* 1 )-g T (v)}X l \ 



f Yi \ Xi {y)dyX, 



\v\<n 
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where \n 1 ' 2 v\ < ix and \n 1 ^ 2 v* l \ < ix when n is large. Both (1A.7j) and (IA.8j) . in con- 
junction with the theorem's assumptions and the finite converging theorem, imply that 



E sup 



v\<J£ 



2.5.1 of 



o(l) for any M > 0. In addition, applying the theorem in Section 



Serflingl ( Il980l ). we have 



1 " 

V^(Qt,Y ~ Qt,y) = fv\Qr,Y) " ~j= ^ MY ~ Q^y) + O p (l) = % 

V n ■ 1 



(1). 



Accordingly, 
1 A 



n 



E 



?T,Y + (Qi-,y— Qt,y) 



fY i \x i (v)dyX i + o p (l) 



i=l ^Qr.y 



1 n 

-(Qr.y - Q T ,y)^ V /y^(Q..y)^ + o P (i; 



- UY f^ T ' Y} A ■ —^MYi- Qr. 



Y) T U p 



oJl) 



Subsequently, using HA.5 j) . flA.6[) . and flA.9j) . we obtain that 
i n 

- J2 MY - Q T , Y )(Xi -X)- qcov T {F, X} 
t=i 
1 n 

-- -=Sy, r {Yi - Q T , Y )(X t - X) - qcov T {E, X}] 

1 n 

= " Qr,v)(X t - Hx\y) - qcov T {F, X}] + 0p (l), 

v 2 = 1 



where //xiy is defined in Subsection 2.2. Since 



\/n(X - /ix) 2 = 
we further have that 



n 

In ^— ' 



1 " 

- <4) = -7= ^2l(X - ^xf - 0- 

v n i= i 



x\ +o P (l). 



(A.9) 



(A.10) 



(A.11) 



Moreover, (lA.lOj) . (1A.11|) . the central limit theorem, and the Cramer- Wold device, lead to 

\ n~ l Eli MY - Qr,y)(A, - X) - qcov T {E, X} ) 
where 



En S 13 

S13 E12 
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and En, E-i?, and Si s are defined in Subsection 2.2. Finally, following the Delta method 
(Ivan der Vaartl . Il998l . Chapter 3), we complete the proof. 



i,z9', 



Proof of Theorem^ We first consider the term <jjq z in qpcof r {y, X|Z}. Let Z* 
X* = Xi — oti — f3[Zi, 61 = and 9\ = (S 1: where (a±, /3[) and (Si , /3() are 

defined in Subsections 2.1 and 2.2, respectively. By the assumptions of this theorem, we 
have that 9 1 = [E{Z*Z*')]- 1 E{Z*X i ), E(Z*X*) = 0, and 

/ 1 n \ ~* 1 n 

yftVi ~ Oi) = ( - E zizr ■ £ z*x; = o p (i). 

According the law of large numbers, we then have that 
1 - 

^iz = ^E(^-^ z *) 2 



i=l 



1 - - /l " 

- Y^&i - + (fi - o x y - J2 

j=l \ i=l 

-m-9 1 y (^E z *x*J 

^(^-^Z^ + OpCn- 1 / 2 ). 



(A.12) 



i=l 



We next consider the numerator in qpcof T {y, X\Z}. For the sake of simplicity, let 
Y* = Yi-a 2 - (3' 2 Zi = Yi- 9 2 Z*, 9 2 = (a 2 ,$)', and 9 2 = (a 2 J 2 )', where Y* is defined 
in the proof of Lemma 2, and (a 2 , f3' 2 ) and (S2, (3' 2 ) are defined in Subsections 2.1 and 2.2, 
respectively. Under the theorem's assu mptions, we emp loy similar techniques to those used 



Koenkerl (120051 ) to show that there exists a unique 



in the proof of Lemma [T] and given in 
9 2 such that E[ip T (Y*)Z*] = and 

1 n 

v^(9 2 - 9 2 ) = {Eify^zmn'w 1 ■ E ^*) z * + ^(i). 

* i=l 

Using a similar method to that for obtaining OA.Sp . we have that 

n ] n 1 n 

- V MY* - 9 2 Z*)X t = - V Vr(l?)* + - E ^ Z - 



i=i 



n 



i=l 



n 



(A.13) 



(A.14) 



where 



g T (Y h Z h 9 2 , 9 2 ) = ip T (Yi - 9 2 Z*) - ip T (Y*) = -[Ifr < 9 2 Z*) - I(Y t < 9 2 Z*)} 
= I[(9 2 - 9 2 )'Z* < Y* < 0] - I[(9 2 - 9 2 )'Z* > Y* > 0]. 



22 



For any v G R q+1 , let 

1 ™ ; f e' 2 z*+n-^ v 'z; 

e„(v) = — y2[g T (Y l ,z l ,e 2 ,e 2 + n-^ 2 v)x l + / fy^xMdyxi 



Applying similar techniques to those for obtaining flA.7[) and (IA.8j) . we can demonstrate 
that 

E[Uv)} 2 = E[g T (Y h Zi, 9 2 , 9 2 + n^V)^] 2 

< < ^^Iv'ZH)} 1 / 2 ■ (M 4 4 ) 1/2 = o(l) 

and, for any 5 > and vi G -R p+1 , 

£ sup \U^i)-U^)\<S-2E[sv^fy i \ Zi ^ i {d , 2 Z* i +y)\X i \]. 

||vi— v||<6 |j/|<T 

This implies that £'supi| v ||< A f |£n(v)| = o(l) for any M > 0. Note that, by (IA.13j) . \fn{6 2 — 
#2) = O p (l). As a result, 



1 n 

1 " r-e> 2 z*+(e 2 -e 2 )>z* 
= — 7^7 Z2 / fv^xMdyXi + p(!) 

1 n 

= -(e 2 - e 2 y ■ -= f YM (o 2 z*)x t z: + o p (l) 

V " • 1 

1 n 

= -E'^ 1 • -= ^r{Y*)Z* + o p (l), 



where E 21 = Eify^xM^XiZ*} and £ 22 = E[f Y . lZi (0 2 Z*)Z*Z*'} are defined in Sub- 
section 2.2. This, together with ( 1A.14|) . results in 

1 n 1 n 

- V MYi - e 2 z;)Xi = - V - 2 z*)pQ - e'^e^z*) + o,^- 1 / 2 ). (A.15) 

rt ^— ' n z — ' 

i=l i=l 



Subsequently, by flA.12j) . (1A.15|) . the central limit theorem, and the Cramer- Wold device, 
we obtain that 



a 2 - a\ 



* |z xlz ]^iV(0,£ 2 ) 
n~ l Er=i MYi ~ d 2 Z*)X t - E[vp T (Y - 6' 2 Z*)X] 



where 

^23 S 25 
E 2 5 E 2 4 
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and S?a, £? 4 , and are denned as in Subsection 2.2. Finally, following the Delta method 
(Ivan der Vaartl . Il998l . Chapter 3), we complete the proof. 



Proof of Theorem^ We first consider the term cr 2 | z in 4>kk,r- Let z\ k _ x = (1, z' t k _i)'- Since 
Ey\ < oo and E[y t — E(y t \ J-^i)] 2 > 0, the matrix ^(z^^z^^) is finite and positive 
definite. Analogous to (IA.12I) . we can show that 



^2 

y| n 

t=k+l 



1 

n. z — ' 

t=k+l 

= E{y t _ k -a x - f3[z^ k _ 1 ) 2 + o p (l). (A.16) 

We next study the numerator of 4>kk,r- Let 8 2 = (0o( r ); <fii(. T ), ■•> 4>p( T ), 0')', and 6> 2 = 
(52, Z?^)'; where is the (k — p) x 1 vector defined in the proof of Lemma 3, and 0.2 and 
02 ar e defined in Subsection 3.1. It is noteworthy that the series {y t } is fitted by model 
(13. ip with order k — 1 and the true parameter vector 6 2 . Accordingly, e tjT = yt — 8' 2 z* tk _ x 
and the parameter estimate of 82 is 82- Then, using ( 1A.19|) in the proof of Theorem HI we 
obtain that 



1 - 

^{82 - 8 2 ) = {E\f t _MKk-i<k -i]}" 1 • - E Met,rKk-i + o v {n' t2 



t=k+l 



Applying a similar approach to that used in obtaining flA.91) . and then using the above 
result, we further have that 



1 n 

~ ^2 ^r{Vt ~ 0'2 z t-k) ~ 4>Ae t , T )]y t _ k 

t=k+l 

= V / ft-i(s)dsy t - k + o p (n 1/2 ) 

t=k+l J[) 

1 " 

= -(82-82)'-- ft-i(0)y t -kzlk-i + o P (n^ 2 ) 

t=k+l 

1 " 

= -Afr^ ■ - M<H,r)*b-i + o P (n- 1/2 ), (A.17) 

t=k+l 

where A\ and E31 are defined as in Subsection 3.1. Subsequently, using similar techniques 
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t=k+l 



to those for obtaining (1A.5|) and the result from equation (1A. 17[) . we obtain that 
1 n 

EiPAvt - «2 - p' 2 z t ,k-i)yt-k 

^ n 1 n 

- E ipr(et,r)yt-k + - E [ipAyt - #2 z tVi) - ip T (e tjT )]y t -k 
n L — ' n z — ' 

t=k+l t=k+l 
1 " 

~ E M^,r)[yt-k - A'^Xfc-i] + o P (n~ 1/2 ). (A.18) 



t=k+l 



Equations (I A. 16[) and (IA. 18[) . together with the central limit theorem for the martingale 
difference sequence, complete the proof of the asymptotic normality of 4>kk,r- From Lemma 
|3l we also have that 4>kk,r = 0. 



Proof of Theorem^ For any v G R p+1 , denote 

n n 

Q(y) = E privt - (^(r) + ™- 1/2 v)'z t y - E - ^)<p) 



t=p+i 



i=P+l 



E Pr( e t,r- n V2v Xp)- E ^Kr)' 
t=p+l t=p+l 



where el T — y t - 
Theorem 3.1 in 



(t>'(T)z* F . Applying fjA.ll) and techniques similar to those in the proof of 



Koenker and Xiaol (120061 ). we can show that 



Q(V) = -V' ■ -= E M<,rKv + E / J ( e *V ^ S ) " ^<r < 0)d8 

* t=p+l t=p+l-' Q 

1 - 1 

= - v '-7^ E ^(e: )T )z* p +-v'£;[/ t _ 1 (0) Z * pZ -]v + Op (l). 
V t= P +i 



Knight 



( 119981 1 . we then have 



Note that Q(v) is a convex function with respect to v. By 
the Bahadur representation as follows, 

1 n 

v^{0(r) - 0(r)} = { J E[/ t _i(0)z*Xp]}~ 1 " -Tjr E <M<r)<p + ^C 1 )- (A. 19) 

V t= P +i 

This, in conjunction with the central limit theorem and the Cramer- Wold device, completes 
the proof. 

Proof of Theorem^ Without loss of generality, we assume that z lp is observable. Then 



e t ,r = Vt~<t> (jK =y t - </>'(t)z* - (0(r) - ^(r))'^ = e* iT - (0(r) - 0(r))'s 
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for 1 < t < n. We first consider the term a 2 in r . By the ergodic theorem and the fact 
that 4>(r) — 4>(r) = O p (n~ 1 / 2 ), we can show that 



and 



t=fc+l t=k+l t=k+l 



n — ' n — ' 

t=k+l t=k+l 
1 n n 

= n E e t - 2(0(r) - (t))' • - £ ^ 

t=fc+l t=fe+l 
1 n 

+ (0(r) - 0(r))' • - £ z*X P ■ (0(r) - 4>{r)) - ft 

t=k+l 

= a 2 e +o p (l), (A.20) 

where a 2 is defined in Subsection 3.2. 

We next consider the numerator of 7>. )T . Using the fact that | ^"=ifc+i ^v(et )T )| < 1, we 
obtain 



1 n 

- ^r(et,r)(et~k,T - fJ>e) 
t=k+l 

1 n ~ ~ 

= - E Vvfo - <Mt)< p )[^ - (0(r) - 4>{t))'zU, p ] + O^n- 1 ) 

t=k+l 
1 n ~i 

= - E ~ ^ ( r Kp) e *-*^ 
f=fc+i 

- (0(r) - 0(r))' ■ - £ Vr(2/t - (r)z* p )z t *_ fcip + o^n" 1 / 2 ). (A.21) 

t=k+l 

Applying similar techniques to those used in obtaining (1A.9j) . we are able to show that 
1 n —/ 

- E ^ T ( yt ~ & ( t Kp) ~ Vv(et, r )]e t _ fe)T 
i 

1 n 

^r(et,r)zt* p + o p (n 1/2 ), 



n 

t=k+l 



n 

t=k+l 



where E 41 is defined in Subsection 3.1 and £51^ = i?[/ t _i(0)e t _fc )r z£' ]. In addition, using 
similar techniques to those in obtaining (1A.5I) and the above result, we further obtain that 



n n 



~ E ^riVt - (P (T)zl p )e t . kiT = - ^ ^r(ei, r ) [et-k,T ~ ^5l,k^41 Z t,p] + ° P ( n V2 )' 



n * — ' ' r n 

t=k+l t=k+l 



2(3 



Analogously, we can verify that 
1 n ~i 



t—k,p 



t=k+l 



The above results, together with (IA.20D . (IA.21D . and the fact that 0(r)-0(r) = O p (n~ 1/2 ), 
imply 



1 1 

TV = , = ■ - V ^r(ei, r )[e 4 _ feiT - Esx^E^z* 



pj + W 



-l/2> 



t=fc+i 



and 



1 1 - 

,x o ' ~ ^r(e t , T )[e t „i, x -S 51 S 41 1 z* p ] + Op(n~ 1/2 ), 



T — T )<7 



t=fc+l 



where e^^x and S 5 i are defined in Subsection 3.2. Subsequently, applying the central limit 
theorem for the martingale difference sequence and the Cramer- Wold device, we complete 
the proof. 
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Table 1: Bias (BIAS), estimated standard deviation (ESD), and asymptotic standard 
deviation (ASD) of the sample quantile correlation qcor r {F, X}. 



BIAS ESD ASD 











h H s 


h B 


3h H s 


0.6/ib 


50 


0.25 


0.0062 


0.1216 


0.1274 


0.1268 


0.1199 


0.1337 




0.50 


-0.0055 


0.1180 


0.1216 


0.1210 


0.1167 


0.1253 




0.75 


0.0023 


0.1221 


0.1286 


0.1280 


0.1199 


0.1364 


100 


0.25 


-0.0045 


0.0828 


0.0867 


0.0860 


0.0836 


0.0891 




0.50 


-0.0032 


0.0792 


0.0835 


0.0829 


0.0816 


0.0847 




0.75 


-0.0029 


0.0824 


0.0873 


0.0863 


0.0836 


0.0897 


200 


0.25 


-0.0012 


0.0598 


0.0601 


0.0596 


0.0586 


0.0607 




0.50 


-0.0000 


0.0585 


0.0580 


0.0577 


0.0571 


0.0582 




0.75 


0.0005 


0.0562 


0.0601 


0.0596 


0.0586 


0.0608 



Table 2: Bias (BIAS), estimated standard deviation (ESD), and asymptotic standard 
deviation (ASD) of the sample quantile partial correlation q^jcor T {Y, X\Z}. 



BIAS ESD ASD 











h H s 


h B 


3h H s 


0.Qh B 


50 


0.25 


-0.0135 


0.1330 


0.1383 


0.1379 


0.1281 


0.1492 




0.50 


-0.0054 


0.1365 


0.1350 


0.1337 


0.1282 


0.1398 




0.75 


0.0029 


0.1378 


0.1407 


0.1401 


0.1299 


0.1487 


100 


0.25 


-0.0094 


0.0901 


0.0972 


0.0963 


0.0922 


0.1004 




0.50 


-0.0026 


0.0931 


0.0943 


0.0935 


0.0912 


0.0959 




0.75 


0.0046 


0.0971 


0.0974 


0.0961 


0.0921 


0.1008 


200 


0.25 


-0.0052 


0.0663 


0.0677 


0.0669 


0.0651 


0.0688 




0.50 


-0.0004 


0.0654 


0.0660 


0.0654 


0.0646 


0.0664 




0.75 


0.0022 


0.0655 


0.0677 


0.0669 


0.0654 


0.0689 
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Table 3: Bias (BIAS), estimated standard deviation (ESD), and asymptotic standard 
deviation (ASD) of the sample QPACF of the observed time series, 4>kk, T , at lags k — 2, 4, 
and 6. 

n t BIAS ESD ASD 



h HS h B 3h HS O.Qh B 



















k = 


2 












50 





,25 


-0 


,0280 


0. 


.1417 


0. 


.1419 


0. 


.1419 


0. 


.1416 


0. 


.1422 







,50 


-0 


,0343 


0. 


1439 


0. 


.1416 





,1416 


0. 


.1416 


0. 


.1417 







.75 


-0 


,0316 


0. 


,1485 


0. 


.1418 





,1419 


0. 


,1416 


0. 


.1421 


100 





.25 


-0 


,0163 


0. 


,1021 


0. 


1001 


0. 


1001 


0. 


1001 


0. 


1002 







.50 


-0 


,0168 


0. 


,1042 


0. 


1001 





,1001 


0. 


1001 


0. 


1001 







.75 


-0 


,0102 


0. 


1009 


0. 


1001 





,1001 


0. 


1001 


0. 


1002 


200 





,25 


-0 


,0107 


0. 


.0732 


0. 


.0707 


0. 


.0707 


0. 


.0707 


0. 


.0707 







,50 


-0 


,0092 


0. 


.0711 


0. 


.0707 





,0707 


0. 


.0707 


0. 


.0707 







,75 


-0 


,0077 


0. 


.0728 


0. 


.0707 





,0707 


0. 


.0707 


0. 


.0707 


















k = 


4 












50 





.25 


-0 


,0344 


0. 


,1434 


0. 


.1438 


0. 


,1439 


0. 


,1421 


0. 


.1457 







.50 


-0 


,0340 


0. 


,1471 


0. 


.1427 





,1424 


0. 


,1427 


0. 


.1437 







.75 


-0 


,0317 


0. 


,1497 


0. 


.1438 





,1439 


0. 


,1421 


0. 


.1456 


100 





.25 


-0 


,0143 


0. 


1032 


0. 


1007 


0. 


1005 


0. 


1002 


0. 


1013 







,50 


-0 


,0172 


0. 


1013 


0. 


1003 





,1002 


0. 


1002 


0. 


1005 







,75 


-0 


,0196 


0. 


1038 


0. 


1007 





,1005 


0. 


1002 


0. 


.1011 


200 





,25 


-0 


,0042 


0. 


0709 


0. 


0709 


0. 


.0708 


0. 


.0708 


0. 


.0710 







,50 


-0 


,0066 


0. 


.0720 


0. 


.0708 





,0708 


0. 


.0708 


0. 


.0708 







.75 


-0 


,0072 


0. 


0703 


0. 


0709 





,0708 


0. 


,0708 


0. 


0710 


















k = 


6 












50 





.25 


-0 


,0278 


0. 


1486 


0. 


1489 





,1489 


0. 


,1483 


0. 


1496 







.50 


-0 


,0356 


0. 


1500 


0. 


.1452 





,1450 


0. 


,1452 


0. 


.1463 







,75 


-0 


,0296 


0. 


.1588 


0. 


.1531 





,1533 


0. 


.1523 


0. 


.1557 


100 





,25 


-0 


,0124 


0. 


.1052 


0. 


.1018 





,1015 


0. 


1004 


0. 


1030 







,50 


-0 


,0197 


0. 


1049 


0. 


1006 





,1004 


0. 


1005 


0. 


.1014 







,75 


-0 


,0189 


0. 


1073 


0. 


.1017 





,1012 


0. 


1004 


0. 


1030 


200 





.25 


-0 


,0103 


0. 


,0741 


0. 


.0712 





,0710 


0. 


,0708 


0. 


.0716 







.50 


-0 


,0112 


0. 


0736 


0. 


0709 





,0708 


0. 


,0708 


0. 


0710 







.75 


-0 


,0105 


0. 


,0727 


0, 


0712 





,0710 


0. 


,0708 


0, 


0715 
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Table 4: Bias (BIAS), estimated standard deviation (ESD), and asymptotic standard 
deviation (ASD) of parameter estimates 0o( r ) and 0i(t). 



n 


T 


Coefficients 


BIAS 


ESD 








ASD 


















h-HS 




h B 


3h HS 


0.6h B 


50 


0.25 


Mr) 


0.0065 





2114 





2285 





2338 





1846 


0.2086 






Mr) 


-0.0487 





1733 





1907 





1943 





1601 


0.1670 




0.50 


M r ) 


-0.0129 





1923 





1957 





2016 





1942 


0.1879 






Mr) 


-0.0471 





1653 





1684 





1736 





1664 


0.1583 




0.75 


M T ) 


-0.0229 





2141 





2287 





2335 





1845 


0.2086 






Mr) 


-0.0495 





1790 





1899 





1938 





1590 


0.1693 


100 


0.25 


M T ) 


0.0016 





1396 





1513 





1553 





1282 


0.1457 






Mr) 


-0.0187 





1208 





1264 





1302 





1102 


0.1189 




0.50 


M r ) 


-0.0039 





1290 





1329 





1370 





1363 


0.1302 






Mr) 


-0.0186 





1151 





1133 





1176 





1169 


0.1092 




0.75 


M r ) 


-0.0118 





1445 





1488 





1534 





1269 


0.1434 






Mr) 


-0.0197 





1227 





1241 





1289 





1091 


0.1164 


200 


0.25 


M T ) 


-0.0010 





0970 





1019 





1052 





0911 


0.1001 






Mr) 


-0.0068 





0839 





0860 





0895 





0785 


0.0831 




0.50 


M T ) 


-0.0046 





0916 





0916 





0943 





1049 


0.0907 






Mr) 


-0.0110 





0765 





0782 





0813 





1082 


0.0769 




0.75 


M r ) 


-0.0099 





0999 





1021 





1055 





0912 


0.1002 






Mr) 


-0.0104 





0840 





0860 





0897 





0787 


0.0833 
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Table 5: Bias (BIAS), estimated standard deviation (ESD) and asymptotic standard de- 
viation (ASD) of the sample QACF of residuals, rjk, T , at k — 1, 3, and 5. 



n QACF BIAS ESD 



ASD 















h H s 


h B 


3h H s 


0.Qh B 
















T = 


0.25 










50 


ri, T 





0195 





0849 





0704 


0.0705 





0699 





0706 




r 3,r 


-0 


0011 





1485 





1478 


0.1479 





1476 





1482 




r 5,r 





0002 





1478 





1497 


0.1496 





1495 





1499 


100 


ri, T 





0075 





0577 





0490 


0.0491 





0489 





0492 




r 3,r 


-0 


0042 





1009 





1008 


0.1008 





1007 





1009 




r 5,r 


-0 


0025 





1041 





1027 


0.1026 





1027 





1029 


200 


n,r 





0035 





0370 





0351 


0.0352 





0351 





0352 




r 3,r 





0008 





0705 





0701 


0.0702 





0701 





0702 




r 5,r 





0005 





0716 





0716 


0.0715 





0716 





0716 
















r = 


0.50 










50 


n,r 





0156 





0779 





0699 


0.0700 





0699 





0704 




r 3,r 


-0 


0199 





1409 





1473 


0.1473 





1473 





1475 




r 5,r 


-0 


0234 





1463 





1492 


0.1491 





1492 





1493 


100 


ri, T 





0094 





0553 





0493 


0.0493 





0493 





0494 




r 3,r 


-0 


0089 





0998 





1008 


0.1007 





1007 





1008 




r 5,r 


-0 


0108 





1022 





1025 


0.1025 





1025 





1026 


200 


r l,T 





0021 





0366 





0351 


0.0351 





0351 





0351 




r 3,r 


-0 


0036 





0693 





0701 


0.0701 





0701 





0701 




r 5,r 


-0 


0049 





0703 
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Table 6: Rejection rate of the test statistic Qbp{K) with K = 6 and the 5% nominal 
significance level. 
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Figure 1: The sample QPACF of the observed time series, 



and 0.8. The dashed lines correspond to ±1.96y f2 3 /n. 
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Figure 2: The time series plot and the sample ACF of the log return (as a percentage) of 
the daily closing price on the Nasdaq Composite from January 1, 2002 to December 31, 
2007. 
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Figure 3: The sample QPACF of daily closing prices on the Nasdaq Composite and the 
sample QACF of residuals from the fitted models for r = 0.2, 0.5, and 0.8. The dashed lines 
in the left and right panels correspond to ±l.9Q\J Vtz/n and ±1.9Qy respectively. 
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